Scalable quantum computing with atomic ensembles 
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Atomic ensembles, comprising clouds of atoms addressed by laser fields, provide an attractive 
system for both the storage of quantum information, and the coherent conversion of quantum infor- 
mation between atomic and optical degrees of freedom. In a landmark paper, Duan et al. (DLCZ) 
[I] showed that atomic ensembles could be used as nodes of a quantum repeater network capable of 
sharing pairwise quantum entanglement between systems separated by arbitrarily large distances. 
In recent years, a number of promising experiments have demonstrated key aspects of this proposal 
[2 03 |U [5] El [7] . Here, we describe a scheme for full scale quantum computing with atomic ensem- 
bles. Our scheme uses similar methods to those already demonstrated experimentally, and yet has 
information processing capabilities far beyond those of a quantum repeater. 



Amongst the more promising schemes for the imple- 
mentation of scalable quantum computing are those in 
which qubits are stored in individual trapped atoms and 
entangled via single photon interference of photons emit- 
ted by the atoms [HI [9]. An appealing aspect of such 
schemes is that once high-fidelity elementary one- and 
two-qubit operations can be demonstrated, it is in prin- 
ciple straightforward to scale to a large number of qubits. 
Notwithstanding recent experimental progress [ID] , this 
approach remains experimentally challenging due to the 
difficulty of trapping and manipulating single atoms, and 
the difficulty of coupling a single atom to a single optical 
cavity mode to improve photon collection efficiency. 

Atomic ensembles provide a promising alternative sys- 
tem. The atomic excitations of a cloud of N identical 
atoms are strongly coupled to the optical field through 
collective enhancement [I], which increases the effective 
atom-photon coupling by a factor of v^TV over the single 
atom case, negating the requirement for a cavity. Fur- 
thermore, it is not necessary for the atoms to be in the 
Lamb-Dicke limit [IT] , significantly reducing the trapping 
and cooling requirements. 

We encode a logical qubit in the collective excitations 
as |0) L = \H) = fft|G) and |1) L ee \V) = W|G). 
Here = N'i^S^ (for £t = #t Q r W) represent 
symmetric collective excitations of the ensemble, where 
= \H(i})(G(i}\ denote the excitation operator for 

atoms, and \G) = |G(i)G(2)...G(jv)) is the state with 
all atoms in the ground state. Two alternative encod- 
ings are possible: an internal-state encoding, where the 
qubit is encoded in two internal atomic levels, in which 
case Hu\/V(i) refer to the two metastable energy levels 
in Fig. alternatively, a dual rail encoding, where the 
qubit is encoded in a pair of identical ensembles, each 
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consisting of atoms with the simpler level structure of 
Fig. [I}}, in which case H/V label the relevant ensemble. 
Note that the two ensembles in Fig.fl}} may simply be two 
spatially distinct regions within the same atomic vapour 
cell, addressed by independent lasers. 

A collective excitation, together with a forward- 
scattered Stokes photon ('excitation' panels of Fig. [I]), 
can be generated by weakly driving one arm of the cor- 
responding Raman transition, yielding the state 

|^> = |G)|vac) + y/pSi S i\G)\vac) + 0(p), (1) 

where |vac) denotes the vacuum state of the optical mode, 
the corresponding photon creation operation operator 
(which could be = h) or v^), and p <C 1 is the prob- 
ability of exciting the ensemble. The third term, corre- 
sponding to multiple excitations, can be made arbitrarily 
small relative to the second term by reducing the excita- 
tion probability p. The collective excitation may be mea- 
sured destructively during 'readout', shown in Fig. [I] by 
driving the reverse transition, resulting in the conversion 
of the atomic excitations into anti-Stokes photons; the 
efficiency of this process can be close to unity owing to 
collective enhancement [7J [IT] . 

Since the qubits are each encoded in separate ensem- 
bles (or ensemble pairs) there is no direct interaction be- 
tween the qubits. Our proposal thus follows DLCZ [1 
and entangles separated ensembles using linear optics 
networks and photodetection on the coherently scattered 
light. Such entangling operations are non-deterministic, 
with success heralded by a particular sequence of pho- 
todetector clicks. To overcome this indeterminism, we use 
a form of measurement based QC [12], in which heralded 
entangling operations can be used to efficiently construct 
an entangled cluster state of many qubits [HI [13] [14] US] • 
Such states are described by a graph comprising a col- 
lection of edges between qubits. The cluster state cor- 
responding to this graph is defined as the state that 
results from initializing each qubit in the state |+) = 
(|0) + |l))/-s/2, and then applying controlled-phase oper- 
ations CZij = diag(l, 1,1,-1) to all qubit-pairs i and j 
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linked by graph edges. Once the cluster state has been 
prepared, universal QC can be implemented by a se- 
quence of single qubit measurements on the state of the 
form sm(0 i)Xi + cos(0 i)Yi. Here Z{ are the Pauli 

operators on the zth qubit, and 6i is a parameter that 
depends on the outcomes of earlier measurements. These 
measurements are implemented by mapping the ensem- 
ble qubit onto a photon ('readout' of Fig. [I]), transform- 
ing the photon polarization using standard linear optical 
elements, then measuring the photon in the h-v basis. 

In the remainder of this paper, we describe a pro- 
cedure for building a cluster state of atomic ensemble 
qubits using linear optical networks and photodetection. 
We first outline a process for fusing two arbitrary cluster 
states together via an optical network that implements a 
heralded controlled-phase operation. We then describe a 
protocol for constructing primitive 3-qubit cluster states. 
Cluster fusion plus 3-qubit cluster states is then sufficient 
to build arbitrary clusters. The protocol is first described 
for the idealized case of perfect 'readout', perfect collec- 
tion and detection efficiency, negligible unwanted excita- 
tions, and negligible decoherence of the atomic ensemble 
qubits. We subsequently argue that the scheme is also 
robust when these idealizations are relaxed. 

For clarity of exposition, we assume that arbitrary lo- 
cal unitary mode transformations can be performed on 
each ensemble (or pair of ensembles) representing a single 
qubit. In fact, such operations may be difficult to imple- 
ment in practice; however, as we describe in the supple- 
mentary information, all such operations can be deferred 
until after the atomic excitations have been mapped onto 
photon states, and thus can be implemented with stan- 
dard linear optical elements. 

We now describe how to apply a controlled-phase (CZ) 
gate between two qubits, labelled 3 and 4 in Fig. [2^i, 
each attached to an arbitrary cluster, and each attached 
to singly linked qubits (labeled 1 and 2). This gate is 
destructive since it consumes qubits 1 and 2, and imple- 
ments a CZ gate between their neighbours. This opera- 
tion serves to forge new links between clusters, thereby 
building larger or more connected graphs. 

The (unnormalized) state in Fig. [2^i may be written 

|$> = (Hi + V}Z 3 )(Hl + T/ 2 t Z 4 )|$')|G) 1 |G) 2 . (2) 

The collective excitations in ensembles 1 and 2 are con- 
verted to photons using the readout pulse of Fig. [I] The 
optical output of 1 passes through a Hadamard gate, H, 
then is mixed with the output of 2 in the optical circuit 
shown in Fig. ^p. Immediately before the photons arrive 
at the photodetectors, the state of the system is given by 

I*) = [(hW a -hlvl)Z A CZ M 

- {h\vl - h\v\)Z^(JL^ (3) 

+ ((hi 2 -hi 2 )(i+z 3 ) 

+ (vt 2 -vt 2 )Z 4 (l-Z s ))]\<t>>)\vzc) ab , 




FIG. 1: Configuration and atomic level structure for (a) 
internal-state encoding: a qubit is encoded in the two internal 
atomic states \H^) and |V(i)), of atoms in a single ensemble 
and (b) dual-rail encoding: a qubit is encoded in the single 
atomic state of atoms in two separate but identical ensem- 
bles, labelled 'if' and 'V. The ensembles may be two dis- 
tinct regions within the same vapour cell, addressed by spa- 
tially separated lasers. When operated in 'excitation' mode, 
a weak, off resonant laser field drives the upward transition 
(straight lines), and a Stokes photon is emitted (wiggly lines). 
In 'readout' mode, the frequencies are reversed, producing an 
anti-Stokes photon in the output mode. 




FIG. 2: a) Protocol for implementing a destructive CZ gate 
between nodes of a cluster state, each designed to have a 
singly-linked node (nodes 1 and 2), with which to build the 
link. Readout pulses are applied to ensembles 1 and 2, convert- 
ing the atomic excitations into photons. The box labelled H is 
a Hadamard gate on the qubit originally encoded in ensemble 
1, which can be implemented via linear optical elements. Note 
that \&) may be a product of two disconnected cluster states, 
so this operation fuses the two together, b) Box B consists of 
a beam splitter (BS) and two polarisation sensitive photode- 
tectors. A successful CZ operation is implemented between 
qubits 3 and 4 when the detected photons are found to have 
opposite polarisations. 
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FIG. 3: Protocol for transforming three two-qubit EME states 
into a three qubit cluster state using polarisation rotators 
and polarising beam splitters (PBSs). Numbered circles rep- 
resent atomic ensembles, double lines indicate EME correla- 
tions, and thin lines represent optical paths. The dashed line 
encloses linear optical elements. Readout pulses are applied to 
qubit s 1, 3 and 5. If each of the three detectors register a single 
photon, ensembles 2, 4 and 6 are projected onto a state that 
is equivalent, up to local transformations, to a three-qubit 
cluster state. 



where a and b label the output modes from the beam 
splitter. The first two terms in represent states in 
which the two photons have different polarisations, and 
the last term represents states in which the two photons 
have identical polarisations. Thus when the photodetec- 
tors click, with 50% chance the photons are found to have 
opposite polarisations, and the state is successfully pro- 
jected onto one that is equivalent (up to local Zi opera- 
tions which must be applied conditional on the outcome) 
to having a CZ gate applied between ensembles 3 and 4, 
resulting in a new cluster state with a link between nodes 
3 and 4, as shown in Fig.[2^i. When photons are registered 
with the same polarisation, the gate fails, and ensemble 3 
is projected onto one of the eigenstates of Z3, which has 
the effect of removing it from the cluster. Thus, in the 
event of a failure, a total of three ensembles are removed 
from the initial state. 

A cluster state can be efficiently constructed with 
this non-deterministic but heralded CZ operation: fail- 
ure of any operation damages the graph, but this can 
be repaired in subsequent steps. Provided an appropri- 
ate strategy is used, the total cost of preparing the state 
is proportional to the size of the cluster [8] [I3j [14] [15] [T6] 

nil m us]. 

With our CZ gate, successful fusion of an n qubit clus- 
ter and an m qubit cluster results in a new cluster of size 
n + m — 2. Thus, in order to grow large cluster states, 
we require an initial supply of 3-qubit cluster states. We 
now describe a two-step recipe for preparing these states, 
starting from an initial supply of ensembles in the \G) 
state. 

The first step is to prepare three copies of an 'effec- 



tive maximally entangled' (EME) state of two ensembles, 
which is given by \EME) id = (H}+V})(V} +H$)\G) [20]. 
These states are 'effectively' entangled, in the sense that 
projecting them onto the subspace with a single excita- 
tion per ensemble results in a maximally entangled state 
of two qubits. EME states of two ensembles (or ensemble 
pairs) may be prepared by first applying a weak excita- 
tion pulse to each ensemble, mixing the forward scattered 
Stokes photons on a 50/50 beamsplitter, and then detect- 
ing the photons with photodetectors. If a single photon 
is registered, the ensembles are left in one of the states 
(Hj ±Hj)\G), depending on which detector clicked. Re- 
peating the procedure with excitation pulses results 
in the state {h\±h]){v} ± V/)|G). By applying appro- 
priate local unitary mode transformations, this state can 
be brought into the form |EME)^ ? j. 

While these states may be useful in certain small scale 
applications [20 , the existence of multiple excitation 
terms are problematic for a scalable quantum computing 
scheme, since they correspond to leakage errors in the 
computation. We overcome this problem with the net- 
work of Fig. [3] This takes as input three EME pairs, and 
conditional on the correct sequence of measurement out- 
comes, outputs a true three-qubit maximally entangled 
cluster state, with only a single excitation per qubit. This 
network was inspired by the observation that a similar 
network, in an all-optical QC context [2T] . also removes 
double excitation terms. Given a supply of ideal EME 
states, the success rate for this step is 1/32 (c.f. [2T]). 

Our destructive CZ gate together with this initial sup- 
ply of 3-qubit clusters is sufficient to build arbitrary clus- 
ter states with a total cost linear in the size of the clus- 
ters. To give an estimate of this cost, our (as yet unopti- 
mized) scheme can produce an TV-ensemble linear cluster 
state with a total of 1536iV/p elementary laser operations 
(see supplementary information). (Note that p should be 
made sufficiently small to reduce the rate of multiple ex- 
citation errors in the state). 

Physical processes which lead to errors in the computa- 
tion include atomic decoherence, losses in linear optical 
elements, imperfect coupling between collective atomic 
and optical modes [7], and imperfect photo-detection 
and dark counts [22 . Fault tolerant quantum compu- 
tation (FTQC) architectures exist for non-determinsitic 
measurement-based schemes [23] [24] [25], so we do not 
address the general question of how to correct all errors 
in our scheme, but note that FTQC can be implemented, 
provided the total error rate lies below the FTQC thresh- 
old, which is around 10 -3 — 10 -2 . 

The dominant sources of error in this proposal - im- 
perfect coupling efficiency, photon loss and detector inef- 
ficiency - are quantified by the effective readout proba- 
bility 77, that an excitation in a collective atomic mode is 
firstly mapped by a 'readout' pulse into the correct op- 
tical mode, then propagates through the optical network 
and is ultimately detected by a photo-detector. A readout 
failure thus leads to a heralded loss error, signified by the 
absence of a photo-detection event. In the supplementary 
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information we draw on previous work [21 to show that 
heralded loss errors can be tolerated provided r] > 2/3, 
which is much less restrictive than FTQC thresholds for 
more general error models. 

This result is particularly promising for our scheme 
since many decoherence processes such as thermal mo- 
tion, atomic dephasing and spontaneous emission simply 
lower the coupling rate between collective atomic and op- 
tical modes [TJ [7] . These errors simply reduce r] (rather 
than producing logical errors), so can be tolerated with 
the very modest threshold for heralded loss errors. This 
robustness is a consequence of the redundant encoding of 
the logical qubits in collective modes of many atoms. 

We have described a scalable scheme to perform quan- 
tum computation with atomic ensembles and linear op- 
tics. The scheme uses similar methods to those used in 



quantum repeater experiments, and yet our proposal has 
information processing capabilities far beyond those of a 
quantum repeater. An important aspect of the scheme 
is the efficient elimination of doubly-excited components 
in the created entangled states. A reasonable near future 
experimental goal is the creation of heralded 3 qubit clus- 
ter states (Fig. [3]). This involves only a moderate increase 
in complexity over existing experiments, and would allow 
interesting applications such as multi-party tests of quan- 
tum mechanics [26], and quantum secret sharing [27] in 
a non-postselected setting. 
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Simon Benjamin, Gerard Milburn, and Terry Rudolph 
for useful conversations. SDB is supported by the EP- 
SRC. PPR thanks the QIPIRC (No. GR/S82176/01) for 
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I. SUPPLEMENTARY MATERIAL: 
OVERCOMING LOSS ERRORS 

In this section we consider in more detail the effect of 
various loss processes within our scheme. In particular, 
the dominant loss processes are due to photodetector in- 
efficiency (that is, the effect of photodetectors failing to 
register a 'click' when a photon enters a detector) and 
the effect of imperfect coupling between the atomic en- 
semble and the correct forward-scattered photon mode. 
Imperfect ensemble-photon coupling arises from a num- 
ber of physical processes. Firstly, there is a fundamental 
limit imposed by the competition between collective cou- 
pling of the ensemble to the forward scattered mode, and 
single-atom spontaneous emission into other free-space 
modes [TJ [7] . Secondly, thermal motion of the atoms can 
reduce the efficiency of the ensemble-photon conversion 
process [7]. In addition, a variety of dephasing and relax- 
ation processes which act on the atoms in the ensembles 
can reduce the effective coupling efficiency [7]. 

In this work, we model detector inefficiency by replac- 
ing each inefficient detector with a perfect one, preceded 
by a beamsplitter with transmissivity rjjj. Similarly, im- 
perfect ensemble-photon coupling processes can be mod- 
eled by assuming the ensemble-photon coupling is per- 
fect, but adding a beamsplitter with transmissivity tje 
on the output of each ensemble. 

An important technique for analysing the effect of 
these losses is to note that, formally, these beamsplitters 
can be commuted through other linear optical elements, 
with the aid of the commutation relations given in [21]. 
To further simplify the analysis we make the assumption 
that each detector has the same efficiency t]d, and that 
the ensemble-photon coupling efficiency takes the same 
value t]e for every ensemble. We also, for the purposes of 
this section, ignore all other imperfections such as detec- 
tor dark counts. This assumption is justified since with 
modern APD detectors, dark count rates are typically 
rather low (~ 50 s _1 ). 

In the remainder of this section, we consider the ef- 
fect of these losses at each stage of our protocol, starting 
at the lowest level (creating two-ensemble EME states), 
and then considering the higher level processes of creat- 
ing 3-qubit GHZ states, and bonding clusters with the 
destructive CZ gate. Our aim is to show that, to a good 
approximation, these errors can lead to an 'independently 
degraded' (ID) error model, where we can form states 
that correspond to initially ideal cluster states which have 
subsequently been subject to uncorrelated qubit loss er- 
rors acting independently on each qubit. The effect of 
such loss errors is to place the affected ensemble qubit 
in the \G) state, regardless of it's initial state. These ID 
states can then be used to perform FTQC with a very 
high threshold [21] , corresponding to a loss probability of 
0.5 per qubit. Throughout this section, we make substan- 
tial use of methods and results due to Varnava et al. [21], 
who discussed the issue of loss tolerance in the context 
of all-optical quantum computing. 
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FIG. 4: Commuting detector loss through the beamsplitter in 
the EME preparation circuit so as to map it to source loss. 
Now the loss in the system can be characterized in terms of 
a single efficiency parameter 77 = t\et\d acting on the source. 



A. Effect of losses in EME state preparation 

Here we consider the effect of loss at the lowest level of 
the protocol, i.e. when we try and make EME states of 
two ensemble qubits. As noted above, there will be two 
types of loss error at this stage of the protocol, charac- 
terized by the detection and coupling efficiencies, 77^ and 
77^, which we can model by beamsplitters with the corre- 
sponding transmissivities (see Fig. |4^l) . The first step in 
the analysis is to commute the beamsplitters tjd to the 
output of each ensemble, leading to an effective model in 
which the detectors are perfect, but each ensemble has a 
beamsplitter of effective transmissivity 77 = tjdVe placed 
at its output (Fig. [4)3). 

Consider a single ensemble, together with the 77- 
beamsplitter at its output. Immediately after an excita- 
tion pulse (acting on, for example, the h Raman transi- 
tion of the ensemble), the combined state of the ensemble 
and the optical fields is [T] 

\^i = T>^( H i h l) n \G)i\0)i, (4) 

n 
n 

(5) 

where we have performed the beamsplitter transforma- 
tion h\ — > y/fjh\ + y/1 — rjhj on the second line, with hi 
denoting the mode reflected from the beamsplitter (i.e. 
the lost light). 

Ultimately, the protocol for creating the EME states 
requires post selection of the case where only one detector 
click is observed (for each round of the protocol) and so 
the relevant terms in the above expression are those with 
only a single photon in the i mode. However, by expand- 
ing out Eq. ([5| it is clear that, in the presence of loss, 
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there will be many such terms, each corresponding to 

a different total number of ensemble excitations, (h}^J . 

Thus, in the final post selected state, we expect to see ad- 
ditional terms which correspond to excess excitations in 
the two ensembles. This can be shown explicitly by con- 
sidering the full network for the EME preparation pro- 
tocol, together with the ^-beamsplitters, tracing out the 
lost modes (i.e. those reflected from the ^-beamsplitters) 
and projecting onto the case where only two detector 
clicks are observed. The result is a state of the form 

PEME,ij = (£ Hi o£ Vi o£ H .o£ v .) ( | EME)ij (EME \ ) , 

(6) 

where the excitation superoperators £s are given, up to 
an overall normalization factor, by 

^HE^^^r^r, (7) 

i>0 

with the excitation operator for the corresponding 
ensemble mode. 

These excitation errors take the ensembles out of the 
logical qubit basis and therefore will lead to errors in the 
computation. Furthermore, it is not obvious that they 
can be detected with certainty in subsequent 'readout' 
stages of the protocol, since by assumption the readout 
and detection efficiencies are less than unity. However, 
by inspecting Eqs. ^ and ^ it is clear that the prob- 
abilities of the error terms are order p(l — rf) or higher. 
Thus the relative contribution of these terms can, in prin- 
ciple, be made arbitrarily small simply by reducing the 
strength of the excitation laser pulse (which controls the 
parameter p). This comes at the expense of decreasing 
the success probability for preparing EME states, which 
incurs a linear increase in the time required to prepare 
these states. Since EME preparation occurs at the lowest 
level of the protocol, and since (given sufficient physical 
resources) many EME states can be prepared in parallel, 
these errors can be strongly suppressed without affecting 
the overall efficiency of the computation. For the remain- 
der of this section we therefore treat the EME states as 
being essentially perfect. In reality, of course, there will 
be some residual multiple excitation errors, but provided 
the magnitude of these is sufficiently small they can be 
dealt with within a more general fault tolerance frame- 
work. 



B. Effect of loss in GHZ state preparation 

We now consider the effect of losses in the GHZ prepa- 
ration step. The inputs to this network are three two- 
qubit EME states, which we assume to be perfect EME 
states, |EME) i} j = (Hj +V/)(V/ +H])\G). Our aim here 
is to demonstrate that, conditional on observing three 
detector clicks in each of the (polarisation resolving) de- 
tectors shown in Fig. |3j the remaining three ensembles 



are projected into an ID-loss state. We make use of sim- 
ilarities between our network for converting EME states 
into GHZ states, and the circuit presented by Varnava et 
al. [21] for generating GHZ states of photonic qubits. 

As in the previous section, we model the photodetec- 
tor and ensemble-photon coupling efficiencies, rjjj and t]e, 
by beamsplitters with the corresponding transmissivities 
placed at the inputs to the detectors and outputs to the 
ensembles, as shown in Fig. [4^i (using the linear optical 
network enclosed by the dashed line in Fig. [3|. The first 
step of the analysis is to make use of the commutation re- 
lations for identical beamsplitters to move the t\e beam- 
splitters forwards through the polarization rotators, and 
the tjd beamsplitters backwards through the network, to 
arrive at the equivalent network in Fig. [4]3, where the 
losses are now represented by the three beamsplitters of 
effective transmissivity r] = tjetJd- 

Varnava et al. showed, via a detailed analysis, that the 
all-optical circuit of Fig. [5^, with imperfect single pho- 
ton sources of efficiency rjs (which is unity for a perfect 
source), leads to a state that is locally equivalent to a 
3-qubit ID-GHZ state on the qubits in modes 2, 3 and 5. 
This state is obtained conditional on a single click being 
observed at each of the detectors. Assuming that the de- 
tectors are perfect, Varnava et al. showed that the output 
state is an ID-GHZ state with local loss rate f = 1 — Vs 
acting on each qubit. 

To make use of this result, first note that the circuit 
considered by Varnava et al. is equivalent to the one 
shown in Fig. (5}d, in which the beamsplitters represent- 
ing source inefficiency are commuted through the first 
row of linear optical elements. Now, if we consider the 
state of the six photons at the position of the broken line 
in Fig. [5)3, it is found that they have the same form as the 
three two-qubit EME states at the input to our GHZ net- 
work, up to an unimportant polarization rotation acting 
independently on each mode. Furthermore, the remain- 
der of the optical network lying below the broken line 
in Fig. is identical to the optical network used in our 
GHZ network, Fig. [3| Thus we can directly apply the re- 
sult of [2T] . which implies that, conditional on a single 
click being observed in each detector in our network, the 
final state of ensembles 2, 4, and 6 is the ID-GHZ state 

pghz,2,4,6 = (£h 2 £v 2 £h 4 £y 4 £h 6 £v 6 ) • • • 

(\GHZ) 2A6 (GHZ\) . (8) 

where the loss superoperators Cs are 

C s (p) = (1 - r)p + r(SpS* + SS*pSS*) . (9) 

The loss rate is given by r — 1 — (2 — tjetjd) -1 - Note 
this is slightly different from the value / determined by 
Varnava et al., owing to the fact that we have not (as yet) 
included the effect of coupling and detector efficiencies for 
the (as yet) unmeasured ensembles 2, 4 and 6. 
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FIG. 5: The equivalence of independent losses (a) before and 
(b) after the first row of PBS's in the GHZ preparation circuit 
of Varnava et al. Single photon sources 1 to 6 emit an H- 
polarised photon. Immediately after the first row of PBSs in 
(b) (at the dark dashed line) the photons are in the state 
|EME)i j2 |EME) 3j 4|EME) 5j 6. 



C. Effect of loss in cluster state preparation and 
measurement based computation 

To complete the demonstration of loss tolerance in our 
scheme, we now consider the effect of loss when building 
large cluster states with our destructive CZ gate, and in 
the single qubit measurement phase of the computation. 

First consider the effect of loss when attempting to 
fuse two ID-cluster states via our destructive CZ gate, 
as shown in Fig. |2j As with the other gates, it is possi- 
ble to model the coupling and detector efficiencies with 
beamsplitters, and again the tjd beamsplitters can be 
commuted back through the optical network such that 
both losses can be represented by a single beamsplitter 
with transmisivity 77 = rjE^lD-, at the output of each en- 



semble. The CZ gate is non-deterministic, with success 
heralded by the observation of two detector clicks in sep- 
arate detectors. Thus, provided the two input states are 
ID-cluster states, the principal effect of the losses is to 
reduce the overall success rate of the gate by a factor 
77/(2 — 77). Note that this rate is a product of the loss 
rate, r, of the input ID-cluster states, together with the 
additional losses incurred in the CZ gate itself. Upon suc- 
cess, the resulting larger cluster state will be an ID state 
with the same loss rate as the input states. This means 
that the CZ gate can be used to build up large ID-cluster 
states of arbitrary shape, and, assuming that additional 
losses incurred while these cluster states are being con- 
structed are negligible, the effective loss rate for these 
clusters will be r = 1 — (2 — 77) _1 . 

The other effect of loss in the CZ gate is that new 'fail- 
ure outcomes' of the gate should now be considered. As 
well as the original failure outcome - observing two pho- 
tons in the same mode - we must also consider the case 
when only one or zero clicks are observed. In this case, 
the input state \&) in Fig. [2] is left in an indeterminate 
(i.e. mixed) state. However, it is generally possible to re- 
cycle significant parts of the input clusters by performing 
a successful Z measurement on qubits '3' and '4'. 

The final stage of our scheme where we must consider 
the effect of loss is the 'readout' phase of the measure- 
ment based computation. In this phase, each qubit in 
the constructed cluster undergoes a single qubit measure- 
ment, which comprises a 'readout' pulse on the ensem- 
ble^), followed by some linear optical operations on the 
output photons, and finally photodetection of the result- 
ing photon state. Clearly, losses will be important here 
since failure to register a photon in one of these mea- 
surements will render the remaining part of the cluster 
in an indeterminate state. Again, we can make use of 
the results of Varnava et al. [21 , who gave an explicit 
protocol for dealing with this loss model. This involves 
encoding each logical qubit in a 'tree structure' compris- 
ing several physical qubits. Provided the total effective 
loss rate for the computation (i.e. the combination of the 
underlying loss rate of the initial cluster, combined with 
the additional losses in the single qubit readout phase) 
is less than 1/2, efficient quantum computation is possi- 
ble. This leads to the final result: losses due to imperfect 
ensemble-photon coupling and detector efficiencies can 
be tolerated provided 



j_ > I 

-7] 2 



(10) 



i.e. provided 77 > 2/3. 



II. SUPPLEMENTARY MATERIAL: 
RESOURCE USAGE 

In order to give an estimate of the resources required 
by our scheme, we calculate the cost of preparation of 
linear clusters in terms of the total number of elementary 
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laser pulses required to build a linear cluster of length N. 
Although linear clusters are not sufficient for performing 
arbitrary computations, they can serve as a resource state 
for building up larger clusters of arbitrary shape, with 
moderate additional overhead. 

There are several stages of the protocol to consider, 
which we consider independently. At the lowest level of 
the protocol we prepare EME states from two separa- 
ble ensembles. EME preparation consists of two rounds, 
where each round is independent and heralded. On aver- 
age this requires 2/p pulses before success. 

Now consider the resource consumption of the proto- 
col for making 3-qubit cluster states of Fig. [3j This net- 
work takes as an input three 2-qubit EME states. It then 
successfully bonds these EME's together into a 3-qubit 
cluster with probability 1/32. Thus the number of op- 
erations (i.e. laser excitation pulses) per 3-qubit cluster 
state is 3 x 32 x 2/p = 192/p. 

Next, suppose our procedure for building up large lin- 
ear clusters is as follows: we have a 'main' cluster, which 
we are attempting to grow, and smaller ancillary clusters 
which we attempt to bond onto the main cluster. Be- 
cause our destructive CZ gate succeeds with probability 
50%, it can be seen that the ancillary clusters must be of 
size at least four qubits if the main clusters is to grow on 
average following each bonding attempt. Thus, we now 
turn our attention to building 4-qubit ancillary clusters. 

Preparing a 4-qubit linear cluster requires two 3-qubit 
linear clusters. Given two such states they can be success- 
fully bonded together to form a 4-cluster with probability 
50%. Thus, the total number of operations per 4-cluster 
is 2 x 2 x 192/p = 768/p. 

Given a resource of ancillary 4-clusters, the 'main' clus- 
ter can be grown at an average rate of 1/2 qubit per 
bonding attempt. So the total number of operations per 
final qubit in the cluster is 2 x TV x 768 /p = 1536N/p. 
Note that, given sufficient experimental resources, many 
of the above steps can be performed in parallel and so 
the actual number of time steps required to grow cluster 
states can be much less than the total number of elemen- 
tary laser operations. 

Here we have adopted an 'incremental' approach to 
building up long clusters, and at several stages we have 
assumed that upon a gate failure, the remaining frag- 
ments of the ancillary clusters are not recycled. Note 
that several alternate strategies could be employed when 
preparing large clusters [I6j [T71 [THJ [19] , which might be 
significantly more resource efficient than the incremental 
approach considered here. 



III. SUPPLEMENTARY MATERIAL: SINGLE 
QUBIT OPERATIONS 



sembles corresponding to a single qubit. Such operations 
may be difficult to implement in practice, particularly 
in the case of dual rail encoding. Here we describe how 
these operations can in fact be deferred until the atomic 
excitations are mapped onto photonic states. 

Such single qubit operations are necessary at several 
stages of the protocol: (1) during the production of the 
EME states; (2) during the production of the three-qubit 
cluster states; (3) after CZ operations between nodes of a 
cluster state; (4) during the 'readout' phase of the mea- 
surement based computation, when each qubit is subject 
to a measurement of the observable sin(^)X^ + cos(^)l^. 

A generic state of n atomic ensembles may be written 

|V> = f(Hl Vl ...Hi V}, ...Hi Vl)\G)® n 

where /(...) is a function of the excitation operators act- 
ing on each ensemble. A local unitary mode transforma- 
tion on the modes of the atomic ensemble, u[ a \ trans- 
forms this state as 

f/f a) = f{H\,vl . . . H'lv'l ...Hi V})\G)® n 

where 

with Uj a 2 x 2 unitary matrix. Subsequently, applying 
the readout operation Ri to the ensembles representing 
the z'th qubit transforms Hj — > h\ and — > v\: 

RiU^Mv^ = f(Hl Vl . . . h'l v'l . . OIG^Ivac); 
with the photon operators given by 

Inspecting the above expressions, it is clear that 
i? i ^ a) |^)|vac) i = /7f ) i? i |^)|vac) i , where U^ o) is a uni- 
tary mode transformation on the zth optical mode with 
the same unitary mode transformation matrix U$ as . 
In other words, local unitary mode transformations of 
atomic ensemble modes can be deferred until after the 
readout pulse, and performed on the optical modes in- 
stead. Such transformations are straightforward to im- 
plement with linear optical elements. The particular se- 
quence of such operations that must be performed gen- 
erally depends on the outcomes of earlier measurements 
in the scheme, and so a modest amount of classical pro- 
cessing and optical switching is also required. 



Throughout the description of the scheme in the main 
section of the paper, we assumed that local unitary mode 
transformations can be performed on each ensemble, or, 
in the case of the dual rail encoding, on each pair of en- 



